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METHOD AND SYSTEM FOR 
IMAGE RECONSTRUCTION 

Background of Invention 

A [0001] This invention relates to imaging methods and systems. More particularly, the 
invention relates to a method and system for image reconstruction. 
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[0002] When a substance such as human tissue is subjected to a uniform magnetic field 
(polarizing field B ), the individual magnetic moments of the spins in the tissue 
attempt to align with this polarizing field, but precess about it in random order at 
| their characteristic Larmor frequency. If the substance, or tissue, is subjected to a 

magnetic field (excitation field B ^ ) which is in the x-y plane and which is near the 
Larmor frequency, the net aligned moment, M , may be rotated, or "tipped", into the 
x-y plane to produce a net transverse magnetic moment M . A signal is emitted by 
the excited spins after the excitation signal B ^ is terminated, this signal may be 
received and processed to form an image. 

[0003] When utilizing these signals to produce images, magnetic field gradients (G G 

x y 

and G ^ ) are employed. Typically, the region to be imaged is scanned by a sequence 
of measurement cycles in which these gradients vary according to the particular 
localization method being used. The resulting set of received NMR signals are 
digitized and processed to reconstruct the image using one of many well known 
reconstruction techniques. 



[0004] 



One method of acquiring an NMR data set from which an image can be 
reconstructed is known as the "spin-warp" technique, which is a variant of the well 
known Fourier transform (FT) imaging technique. The spin-warp technique is 
discussed in an article entitled "Spin-Warp NMR Imaging and Applications to Human 
Whole-Body Imaging" by W.A. Edelstein et al., Physics in Medicine and Biology, Vol. 
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25, pp. 751-756 (1980). It employs a variable amplitude phase encoding magnetic 
field gradient pulse prior to the acquisition of NMR spin-echo signals to phase encode 
spatial information in the direction of this gradient. In a two-dimensional 
implementation (2DFT), for example, spatial information is encoded in one direction 
by applying a phase encoding gradient (G ) along that direction, and then a spin- 
echo signal is acquired in the presence of a readout magnetic field gradient (G ) in a 

x 

direction orthogonal to the phase encoding direction. The readout gradient present 
during the spin-echo acquisition encodes spatial information in the orthogonal 
direction. In a typical 2DFT image acquisition, a series of pulse sequences is 
performed in which the magnitude of the phase encoding gradient pulse G in the 

y 

pulse sequence is incremented ( A G ). The resulting series of views that is acquired 
during the scan form an NMR image data set from which an image can be 
reconstructed. The acquisition of each phase encoded view requires a finite amount of 
time, and the more views that are required to obtain an image of the prescribed field 
of view (FOV) and spatial resolution, the longer the total scan time. 

[0005] Many technical developments in the field of MR imaging aim to reduce data 

acquisition time. One such technical development is known as parallel imaging, in 
which multiple receive coils are used to reduce scan time. 

[0006] In parallel imaging techniques, the reduction in scan time comes from omitting k- 
space lines or, equivalently, spacing the lines further apart. However, the reduction in 
k-space lines has the effect of reducing the FOV, thus producing image wrapping or 
aliasing if the object extends outside the reduced FOV in the phase encoding 
direction. 

[0007] 

The aliasing associated with parallel imaging may be eliminated using one of two 
general techniques. In one technique, the missing k-space lines are synthesized using 
the receive coil sensitivity (B1 field) spatial information, followed by a Fourier 
transform to reconstruct an image without aliasing. This technique is called SiMultan 
eous Acquisition of Spatial Harmonics (SMASH), and is described by Sodickson, et al., 
"Simultaneous Acquisition Of Spatial harmonics (SMASH)", Magnetic Resonance In 
Medicine 1997,38:591-603. In a second technique, the aliased images reconstructed 
from each coil are combined using the receive coil sensitivities to cancel the aliasing. 
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This technique is called sensitivity encoding (SENSE) and is described by K.P. 
Pruessmann, et al., "SENSE: Sensitivity Encoding for Fast MRI", Magnetic Resonance in 
Medicine 42, 952-962 (1 999). Although these are the two most common parallel 
imaging methods, other generalized parallel imaging methods combine aspects of 
both techniques. 

[0008] In the SENSE technique, the idea is to reduce acquisition time by increasing the 
step size (AG ) between phase encoding views, or equivalently, by reducing the 
FOV. In either case, the total number of views is reduced with a consequent reduction 
in scan time. If the object extends outside the reduced FOV, however, aliasing or 
wrap-around occurs in the phase encoding direction. The SENSE technique removes 
this aliasing by using knowledge of the surface coil receive field (also called 
sensitivities) to find the unaliased spin distribution. This knowledge is obtained by 
71 performing a sensitivity calibration scan prior to the reduced FOV (SENSE) scan. For 

example, the sensitivity calibration data may be obtained using combined signals 
from the surface coils as described in co-pending U.S. Pat. Appln. Serial No. 
09/851 ,775 filed on May 9, 2001 and entitled "Calibration Method For Use With 
53 Sensitivity Encoding MRI Acquisition". 
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[0009] During both the SENSE scan and the calibration scan, the patient may be required 
to hold his or her breath to minimize respiratory motion artifacts. If the patient and/or 
coil move in the phase encoding direction between the calibration and SENSE scans, 
for example if the patient were to intake a different volume of air during each of the 
scans, artifacts due to uncorrected aliasing can result. 



[0010] 



One method of reducing artifacts due to patient motion between the calibration 
and SENSE scans is to extend the sensitivity measurement outside the boundaries of 
the object boundaries detected in the calibration scan using, for example, linear 
extrapolation. In conjunction with extending the sensitivity, the object boundaries are 
also extended to be outside those measured in the calibration scan. If the sensitivity 
and edges are extended sufficiently far, the sensitivity will be estimated at all pixels 
within the object at its position during the SENSE scan, and in addition, the number of 
aliased replicates will be greater than or equal to the actual number of replicates at all 
pixels. For typical scan prescriptions and typical respiratory motion, it has been found 
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that edge and sensitivity extension by 1 0 to 20 pixels is sufficient. This enables 
accurate aliasing unwrapping at all pixels, even if the object moved after the 
calibration scan. However, the final image noise will be higher at pixels where the 
number of aliased replicates is overestimated. The increased noise frequently covers a 
large enough area to be noticeable and objectionable. 

[001 1] If the number of aliased replicates is overestimated, an adaptive, pixel-dependent 
regularization of the SENSE reconstruction has been used to reduce image noise. 
Adaptive regularization is a well-known numerical method used to reduce noise in the 
solution of inverse problems. However, one drawback to adaptive regularization is 
that it requires much iteration to determine the optimum regularization and therefore 
requires a long reconstruction time. 



yj [001 2] Another method of reducing artifacts due to patient motion between the 

calibration and SENSE scans is to combine both scans while the patient holds a single 
breath. Yet another method of reducing artifacts due to patient motion between the 
calibration and SENSE scans is to allow the patient to breathe freely while taking 



Ml 
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multiple averages during the calibration scan. These methods, however, increase the 



fU 



SENSE scan time, thus partially negating some of the advantages of the SENSE 
technique. 

Summary of Invention 

[001 3] The above-described drawbacks and deficiencies are overcome or alleviated by a 
method for reconstructing image of an object, the method comprising: acquiring a 
calibration image of the object; finding an edge of the object in the calibration image; 
calculating a first number of aliasing replicates at a pixel location in the reduced field 
of view image using the edge found in the calibration image; unwrapping the reduced 
field of view image using the first number of aliasing replicates to create a first 
unwrapped image; finding the edge of the object in the first unwrapped image; 
calculating a second number of aliasing replicates at the pixel location in the reduced 
field of view image using the edge found in the first unwrapped image; and 
unwrapping the reduced field of view image using the second number of aliasing 
replicates to create a second unwrapped image. 
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[001 4] The above discussed and other features and advantages of the present invention 
will be appreciated and understood by those skilled in the art from the following 
detailed description and drawings. 

Brief Description of Drawings 

[001 5] Referring to the exemplary drawings wherein like elements are numbered alike in 
the several Figures: 

[001 6] Fig. 1 is a block diagram of an MRI system which employs the present invention; 

[001 7] Fig. 2 is a schematic representation of a patient in the MRI system of Fig. 1 
surrounded by an array of four local coils; 
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[001 8] Fig. 3 is a flow chart illustrating a method of image reconstruction; 
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H [001 9] Fig. 4 is a pictorial representation of an object scanned with the MRI system of Fig. 
1 , illustrating phase encode columns of a calibration image; 



[0020] Fig. 5 is a pictorial representation of a of an object scanned with the MRI system 

•41 

£l3 of Fi 9- 1 » illustrating a phase encode column of a calibration image; 

l| 

Q [0021] Fig. 6 is a pictorial representation of a of an object scanned with the MRI system 
of Fig. 1 , illustrating a phase encode column of a SENSE image; 

[0022] Fig. 7 is a flow chart illustrating an alternative embodiment of the method of 
image reconstruction; 

[0023] Fig. 8 is a pictorial representation of an object that is bigger than the normal FOV; 

[0024] Fig. 9 is a pictorial representation of the object of Fig. 8 scanned with the MRI 
system of Fig. 1 , illustrating a phase encode column of a SENSE image; and 

[0025] Fig. 1 0 is a pictorial representation of a combination of a calibration image and a 
SENSE image. 

Detailed Description 

[0026] 

Referring first to Fig. 1 , there is shown an embodiment of a MRI system which 
incorporates the present invention. The operation of the system is controlled from an 
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operator console 100 which includes a keyboard and control panel 102 and a display 
1 04. The console 1 00 communicates through a link 1 1 6 with a separate computer 
system 107 that enables an operator to control the production and display of images 
on the screen 104. The computer system 107 includes a number of modules that 
communicate with each other through a backplane. These include an image processor 
module 106, a CPU module 108 and a memory module 113, known in the art as a 
frame buffer for storing image data arrays. The computer system 1 07 is linked to a 
disk storage 1 1 1 and a tape drive 1 1 2 for storage of image data and programs, and it 
communicates with a separate system control 1 22 through a high speed serial link 
115. 
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[0027] The system control 1 22 includes a set of modules connected together by a 

backplane. These include a CPU module 11 9 and a pulse generator module 121 which 
connects to the operator console 1 00 through a serial link 1 25. It is through this link 
125 that the system control 122 receives commands from the operator which indicate 
the scan sequence that is to be performed. The pulse generator module 121 operates 
the system components to carry out the desired scan sequence. It produces data 
which indicates the timing, strength and shape of the RF pulses which are to be 
produced, and the timing of and length of the data acquisition window. The pulse 
y generator module 121 connects to a set of gradient amplifiers 1 27, to indicate the 

timing and shape of the gradient pulses to be produced during the scan. The pulse 
generator module 121 also receives patient data from a physiological acquisition 
controller 1 29 that receives signals from a number of different sensors connected to 
the patient, such as ECG signals from electrodes or respiratory signals from a bellows. 
And finally, the pulse generator module 121 connects to a scan room interface circuit 
1 33 which receives signals from various sensors associated with the condition of the 
patient and the magnet system. It is also through the scan room interface circuit 1 33 
that a patient positioning system 1 34 receives commands to move the patient to the 
desired position for the scan. 



[0028] 



The gradient waveforms produced by the pulse generator module 1 21 are applied 

to a gradient amplifier system 1 27 comprised of C , G and G amplifiers. Each 

x y z 

gradient amplifier excites a corresponding gradient coil in an assembly generally 
designated 139 to produce the magnetic field gradients used for position encoding 
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acquired signals. The gradient coil assembly 1 39 forms part of a magnet assembly 
141 which includes a polarizing magnet 140 and a whole-body RF coil 1 52. 



[0029] A transceiver module 1 50 in the system control 1 22 produces pulses which are 

amplified by an RF amplifier 151 and coupled to the RF coil 1 52. The resulting signals 
radiated by the excited spins in the patient may be sensed by the same RF coil 1 52, 
but in the preferred embodiment a local coil array indicated at 1 54 is employed. The 
local coil array 1 54 is positioned around the patient anatomy to be imaged and it 
includes four separate coil elements (shown in Fig. 2) which receive the NMR signals 
that are produced during the scan. These NMR signals are separately amplified by 
pre-amplifiers 1 53 and applied to the inputs of four separate receivers in the 
transceiver module 1 50. The amplified NMR signals are demodulated, filtered, and 
digitized in the receiver sections of the transceiver 1 50 to produce four separate k- 
space data sets. 



M= [0030] The local coil array 1 54 may be similar to that disclosed in U.S. patent No. 

:;; 

ui 



4,825,162. Referring to Fig. 2, each multi-coil set consists of two flexible paddles 10 
and 1 2, each with two surface coil elements fabricated with etched copper on a 
flexible plastic substrate, for a total of four coil elements. Each coil element is a 
rectangular loop, approximately 12x25 cm and adjacent elements are overlapped to 
minimize mutual inductance. The elements are tuned using capacitors distributed 
around each element, and interface circuitry is added to permit detuning of the 
element during RF transmission. The coils are positioned anterior and posterior to the 
imaging volume as depicted in Fig. 2. It should be apparent to those skilled in the art 
that the present invention is not limited to any particular local coil array and that 
many alternative local coils are commercially available and suitable for this purpose. 

The NMR signals picked up by the coil array 1 54 are digitized by the transceiver 
module 1 50 and transferred to a memory module 1 60 in the system control 1 22. 
When the image scan is completed the four resulting k-space data sets are processed 
as will be described in more detail below. This data is conveyed through the serial link 
1 1 5 to the computer system 1 07 where it is stored in the disk memory 1 1 1 . In 
response to commands received from the operator console 1 00, this data may be 
archived on the tape drive 1 1 2, or it may be further processed by the image processor 
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106 and conveyed to the operator console 100 and presented on the display 1 04. 
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[0032] The above-described MRI system may be used to perform many prescribed scans 
on a patient positioned in the bore of the magnet 141 . Such prescriptions identify the 
type of imaging pulse sequence to be used (e.g., GRE, FSE, EPI, etc.) and the particular 
scan parameters for the selected pulse sequence (e.g. TR, TE, flip-angle). The present 
invention may be used with any such prescribed scan. 

[0033] Referring particularly to Fig. 3, a method of image reconstruction, which may be 
employed by the system of Fig. 1 , is shown generally at 200. The scan parameters 
which define the imaging volume are input by the operator along with other scan 
parameters. The scan parameters are employed to determine the parameters of a 
calibration pre-scan as indicated at block 202. In Fig. 4, three such prescribed phase 
encode columns are shown at 1 4. As represented in Fig. 5, the FOV of these phase 
encode columns is larger than the object 16 being imaged. In the preferred 

tj 

embodiment the calibration pre-scan employs the surface coil elements 10 and 12 of 
^ Fig. 2 to acquire data from which their sensitivity is determined as disclosed in the 

ff above-cited co-pending U.S. Pat. Appln. Serial No. 09/851,775. Although many 

y different imaging pulse sequences may be used to perform the calibration pre-scan, a 

fast, 2D gradient-recalled echo pulse sequence may be employed because of its speed 
and other well known qualities. A flip-angle of 45 ° , a TR of 1 00 msec and a TE of 6 
msec is prescribed for the calibration pre-scan. The entire imaging volume is covered 
using contiguous slices from 5 to 1 0 mm in thickness and with 1 28 x 128 voxels in 
each slice. These parameters are lower in resolution than the typical imaging 
prescription, but the resolution is sufficient to produce a quality sensitivity matrix S 
while limiting the calibration pre-scan time. 

[0034] In block 204, a calibration scan is performed to obtain calibration data for each 
coil. This is done with the patient positioned in the MRI system. After the last 
calibration slice has been acquired, the prescribed image (SENSE) scan is performed, 
as indicated in block 206. As indicated above, any imaging pulse sequence may be 
used, but RF excitation is performed using the RF body coil 1 52 and NMR signal 
acquisition is performed using the local coil array 1 54. 

[0035] The SENSE technique allows for a faster scan by reducing the field of view (FOV) 
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along one or more phase encoding axes by a reduction factor R. If the field of view is 
reduced in a phase encoding axis by a factor of R, the scan time is also 
correspondingly reduced by a factor of R. However, now the reconstructed image is 
aliased or replicated in the phase encoding direction at multiples of FOV/R=D and 
aliasing replicates now overlap with resulting loss of diagnostic utility. For example, 
Fig. 6 shows a single phase encode column of the SENSE image with a reduced FOV, 
indicated as D. The number of aliasing replicates for each pixel within various regions 
of the image is shown. For the SENSE image to be of diagnostic utility, the image must 
be unwrapped (un-aliased). In Fig. 6, and in the following discussion, the phase 
encoding gradient axis is indicated as the / -axis, however, it should be understood, 
that any phase encoding axis is possible and the invention is not limited to a specific 
gradient axis direction. 

[0036] Referring again to Fig. 3, after the SENSE data has been acquired, the next step in 
method 200 is to reconstruct each aliased SENSE image *,(y) produced by each local 
coily=0, 1, 2, ... N-l, as indicated in block 208. This may include performing a 2D 
Fourier transformation on the aliased SENSE image data and organizing the complex 
intensity values for each location xinto an array I where y \s the phase encode 
location at the given x: 



0) 
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[0037] The calibration data that acquired in block 204 is used to reconstruct a set of 
calibration images *f(y) for each surface coily-0, 2, ... N-l in the phase encoding 
gradient axis y, as indicated at block 21 0. This reconstruction may be a two- 
dimensional, complex Fourier transformation which yields a 1 28 x 1 28 array of 
complex pixel intensity values for each slice acquired by each local coil. 



[0038] 



The sensitivity S . (y) for each local coil j is then calculated at block 21 2. While coil 
J 

sensitivity may be calculated using any known technique, two exemplary 
embodiments are as follows. In one embodiment the local coil sensitivity S . (y) is 
calculated as follows: 
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(2) where N is the number of local coils and rf{y) is a complex image acquired over 
the full FOV of the imaging volume with surface coil j. 
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[0039] In a second embodiment, the sensitivity S . (y) is defined as: 



I Too 

(3) The sensitivity definition in equation 3 has the advantage that the complex 
magnetization completely drops out of the definition, thereby resulting in more 
accurate phase information in the reconstructed image. 



[0040] Both of these methods for calculating the local coil sensitivity S . (y) in essence 
UJ j 

calculate the ratio of the calibration image intensity of the local coil itself and the sum 



of all the local coil calibration image intensities. 



[0041] The calibration images are not necessarily aligned with their corresponding 
SENSE image plane because the different images may be acquired with different slice 
orientations and different image resolution. The next step, therefore, is to align the 
calibration object images if(y) with their corresponding SENSE image planes, as 
indicated at block 214. This can be accomplished by linear interpolation of calibration 
object image pixels. 

[0042] The object edges, which are indicated in Fig. 5, are then located in the calibration 
image, as indicated at block 21 6 of Fig. 3. This may be performed by searching inward 
from the edges of the calibration image along the phase encoding direction until a 
pixel having an intensity above the established threshold is located. All pixels 
encountered before this location are set to zero, and if the center of the image is 
reached, the search is stopped along that path. As a result, a calibration image is 
produced in which all pixels outside the object 16 are set to zero. 

[0043] 

In block 21 8, to help compensate for possible movement of object 16 between the 
calibration and SENSE scans, the sensitivity measurement may be extended outside 
the object boundaries using, for example, linear extrapolation. In block 217, the 
object boundaries may also be extended outside those measured. If the sensitivity 
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edges are extended sufficiently far, the sensitivity will be estimated at all pixels within 
the object at its position during the SENSE scan. In addition, calculation of the number 
of aliased replicates will be greater than or equal to the actual number of replicates at 
all pixels. For typical scan prescriptions and typical respiratory motion, it has been 
found that edge and sensitivity extension by 10 to 20 pixels is sufficient to enable 
accurate aliasing unwrapping at all pixels, even if the object moved between the 
calibration and SENSE scans. 

[0044] The next step indicated at block 220 is to calculate the number of aliased 

replicates, A, using the FOV reduction factor R and the edge locations. Once the object 
edges are determined for each phase encoding column, the number of overlapped 
aliasing replicates for each pixel in the column is determined as follows. Let the pixel 
locations of the upper and lower edges be e y and e respectively. Translate the left 
edge to the right in increments of the alias distance D until the right edge of the 
yj unreduced field of view is crossed. For each pixel count the number of locations, 

including the original untranslated location, for which the pixel is above the edge. Let 
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that number be n Translate the right edge to the left in increments of the alias 

i. n L. 



ru 



distance D until the left edge of the unreduced FOV is crossed. For each pixel y, count 
the number of locations, including the original untranslated location, for which the 
pixel is below the edge. Let that number be n The number of aliased replicates A at 
pixel y\s n ^ + n ^ - /. 

[0045] Method 200 then continues to block 221, where the sensitivity matrix S, is 

populated with separate coil sensitivities S ^ (y), the number of aliased replicates A, 
and the reduced FOV D: 



S = 



(4) 



So(y) S 0 (y+D)„.S 0 (y+(A-l)D) 
S|(y) S l (y+D)...S I (y+(A-l)D) 

S N _i(y) S N . l (y+D)...S N . ! (y+(A-l)D)^ 



[0046] ln block 222, the SENSE images ',00 are then unwrapped to create an unwrapped 
SENSE image (proton distribution image) M. This may include applying the equation: 
M={(s*sy l s*]i 

(5) where S is the sensitivity matrix from block 220, I is the aliased SENSE image array 
of block 208, and S is the complex conjugate transpose of the sensitivity matrix S. 
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[0047] At process block 224 it is determined whether the current iteration is the last 

iteration to be performed. More than two iterations may be performed; however, it is 
anticipated that only two iterations are required for optimal image quality. The 
number of iterations may be predetermined. Alternatively, the number of iterations 

may be determined by comparing the image from the current iteration with the image 

-? 

from the previous iteration. When a suitable measurement of the difference is 
sufficiently small (such a measurement could be the sum of squares of the differences 
for all pixels in the two images), the iteration may be stopped . If the present iteration 
is the last iteration, method 200 ends, and the proton distribution image M is applied 
in a manner known in the art. 



O [0048] If, at block 224, the present iteration is not the last iteration, then method 200 
m continues to block 226, where the object boundary (edge) is located in the proton 

Wj distribution image M from block 222. This may be performed by searching inward 

Ly from the edges of the image along the phase encoding direction until a pixel having 

an intensity above the established threshold is located. 
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[0049] In one embodiment, after the edges have been located, method 200 then returns 
to block 220, where the number of aliased replicates, A, is calculated using the FOV 
reduction factor R and the new edge locations determined in block 226. The method 
used for determining A may be the same as that described for use in block 220. 

[0050] In an alternative embodiment, which is indicated in phantom on Fig. 3, method 

proceeds from block 226 to block 219, where the edge locations determined in block 
226 are extended prior to calculating the number of aliased replicates, A. Extending 
the edge location helps to ensure that any errors in edge detection do not result in 
underestimation of A. For example, the edge locations may be extended by one or two 
pixels, which would increase image noise only in a small area of the image. 

[0051] In either embodiment, after A has been calculated, method 200 proceeds to block 
221 where A is input into the sensitivity matrix S, along with the separate coil 
sensitivities S (y) determined in block 212, the reduced phase encoding FOV, D. In 
block 222, the sensitivity matrix S from block 221 is applied to unwrap the SENSE 
images ',00 from block 208 using, for example, equation (5). Method 200 continues 
until the last iteration is performed. 
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[0052] Referring to Fig. 7, an alternative embodiment of method 200 is shown. The 

embodiment of Fig. 7 is substantially similar to the embodiment of Fig. 3, with a block 
230 added between previously-described blocks 224 and 226. In the embodiment of 
Fig. 7, method 200 is adapted for scan prescriptions where the object is not 
completely contained within the FOV after the image is unwrapped. This would occur, 
for example, in scan prescriptions where the FOV of the normal (non-calibration) 
image is smaller than the object being scanned, as shown in phase encode column of 
Fig. 8. The resulting SENSE image would have a reduced FOV, as indicated as D in the 
phase encode column of Fig. 9, and would yield overlapping aliasing replicates. 
Unwrapping of the aliased SENSE image would yield an image in which the object 
edges are not included, and, therefore, the unwrapped image would not be of value in 
the edge detection block 226. 
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3 [0053] In block 230, the first reconstructed SENSE image is superimposed on the 

calibration image from block 210 that is used to reconstruct the first SENSE image. As 
shown in the phase encode column depicted in Fig. 10, the calibration scan uses a 
large FOV that will normally contain the entire object and, by combining the SENSE 
and calibration images, the edges can be determined. The SENSE and calibration 
images may have different scan parameters, leading to small intensity discontinuities, 
which are exaggerated in Fig. 10. This combination of the SENSE and calibration 
images from block 230 is then used in block 226 to find the edges. Method 200 
Continues as described with reference to Fig. 3. 

[0054] The iterative reconstruction method 200 described herein reduces artifacts due to 
uncorrected aliasing that can result from patient motion between the calibration and 
SENSE scans. In method 200 edge boundaries are estimated using an unwrapped 
SENSE image. As a result, the edge boundaries are correctly estimated even where the 
object moves between the calibration and SENSE images. Because the edge locations 
are accurate, calculation of the number of aliased replicates at each pixel is also 
accurate, resulting in full unwrapping of all aliasing together with the lowest possible 
noise at each pixel. Method 200 accomplishes this without the long scan time or 
reconstruction time associated with known methods. 



[0055] 



The present invention can be embodied in the form of computer-implemented 
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processes and apparatuses for practicing those processes. The present invention can 
also be embodied in the form of computer program code containing instructions 
embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any 
other computer-readable storage medium, wherein, when the computer program code 
is loaded into and executed by a computer, the computer becomes an apparatus for 
practicing the invention. The present invention can also be embodied in the form of 
computer program code, for example, whether stored in a storage medium, loaded 
into and/or executed by a computer, or transmitted over some transmission medium, 
such as over electrical wiring or cabling, through fiber optics, or via electromagnetic 
radiation, wherein, when the computer program code is loaded into and executed by a 
computer, the computer becomes an apparatus for practicing the invention. When 
implemented on a general-purpose microprocessor, the computer program code 
segments configure the microprocessor to create specific logic circuits. 

[0056] While the invention has been described with reference to a preferred embodiment, 
it will be understood by those skilled in the art that various changes may be made and 



w equivalents may be substituted for elements thereof without departing from the scope 

i'J. jj 

Q of the invention. In addition, many modifications may be made to adapt a particular 

m 

;^ situation or material to the teachings of the invention without departing from the 

til essential scope thereof. Therefore, it is intended that the invention not be limited to 

the particular embodiment disclosed as the best mode contemplated for carrying out 
this invention, but that the invention will include all embodiments falling within the 
scope of the appended claims. Moreover, the use of the terms first, second, etc. do 
not denote any order or importance, but rather the terms first, second, etc. are used 
to distinguish one element from another. 
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